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Summary. The propagation of cellularly stabilized thermonuclear flames is inves- 
tigated by means of numerical simulations. In Type la supernova explosions the 
corresponding burning regime establishes at scales below the Gibson length. The 
cellular flame stabilization — which is a result of an interplay between the Landau- 
Darrieus instability and a nonlinear stabilization mechanism — is studied for the 
case of propagation into quiescent fuel as well as interaction with vortical fuel flows. 
Our simulations indicate that in thermonuclear supernova explosions stable cellu- 
lar flames develop around the Gibson scale and that a deflagration-to-detonation 
transition is unlikely to be triggered from flame evolution effects here. 



1 Introduction 

The standard model of Type la supernovae (SNe la) describes these astro- 
physical events as thermonuclear explosions of white dwarf (WD) stars [1]. 
In our study, we refer to the specific scenario (for a review on SN la explo- 
sion models see [2]), where the WD consists of carbon and oxygen and the 
thermonuclear reaction propagates in form of a flame that starts out in the 
so-called deflagration mode. That is, the flame is mediated by microphysical 
transport processes and its burning velocity is subsonic. One key ingredient 
in such a SN la model is the determination of the effective propagation veloc- 
ity of the deflagration flame. The so-called laminar burning speed s;, i.e. the 
propagation velocity of a planar deflagration flame, is much too low (a few 
percent of the sound speed in the unburnt material) to explain powerful SN 
la explosions. The solution to this problem is provided by the concept of 
turbulent combustion. Instabilities on large scales result in the formation of 
a turbulent cascade, where large-scale eddies decay into smaller ones thereby 
transporting kinetic energy from large to small scales. Interaction of the flame 
with those eddies wrinkles the flame front and enlarges its surface. This is 
equivalent to an increase in the net fuel consumption rate and hence causes 
an acceleration of the flame. 

Recent SN la explosion models on scales of the WD could show that in 
this way enough energy can be released to gravitationally unbind the star 
[3, 4]. However, these models cannot resolve all relevant length scales down 
to the flame width and therefore have to rely on assumptions on the physics 
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on unresolved small scales. The goal of our small-scale simulations of the 
propagation of the thermonuclear flame is to test those assumptions and 
eventually to reveal new physics that additionally needs to be included in the 
SN la models. The most significant assumptions of the large scale models by 
Reinecke et al. (e.g. [3]) are that the energy input into the turbulent cascade 
originate solely from large scales and the flame be stable on unresolved scales. 

A fundamental feature of the turbulent cascade is that in a certain in- 
termediate scale range, the turbulent velocity fluctuations decrease mono- 
tonically with smaller length scales of the eddies. This is captured by the 
corresponding scaling law. From that effect it is obvious, that there must 
exist a cutoff scale, below which flame propagation is not affected by the tur- 
bulent cascade anymore. Below this Gibson scale, the velocity fluctuations 
are so small, that the flame burns faster through the eddies than they can 
deform it. Thus the flame here propagates through "frozen turbulence" . Our 
simulations aim on the study of effects around and below the Gibson scale. 
Except for the very late stages of the SN la explosion, the Gibson scale is well- 
separated from the width of the flame and the flame may safely be regarded 
as a discontinuity between burnt and unburnt material. 

2 Theoretical background 

What are the effects that dominate flame propagation below the Gibson scale? 
Here, two counteracting effects become important. The first is the Landau- 
Darrieus (LD) instability [5, 6]. It is based on a hydrodynamical effect and 
would completely prevent flame stability on the scales under consideration. 
This would have drastic impact on SN la models. However, there exists a 
second effect of purely geometrical origin [7], which balances the LD insta- 
bility. Once the perturbations have grown to a critical size, former recesses 
of the flame front will develop into cusps, which possess propagation veloci- 
ties slightly exceeding the burning speed of the rest of the front. This effect 
gives rise to a stable cellular shape of the flame. Therefore, the regime of 
flame propagation that will be studied in the present contribution is termed 
cellular burning regime. 

In connection with SN la explosions, this burning regime was studied by 
[8, 9] and [10] demonstrated its existence by means of a full hydrodynamical 
simulation. The questions that arise here are: Can the cellular flame stabiliza- 
tion break down under certain conditions (e.g. low fuel densities or interaction 
with turbulent velocity fields)? If so, could the flame itself generate turbu- 
lence on small scales (in contradiction to the assumption of large-scale SN 
la models) and thereby actively accelerate? This conjecture was put forward 
under the name active turbulent combustion by [11] on the basis of [12]. The 
questions are closely connected to the search for a mechanism providing a 
deflagration-to-detonation transition (DDT) of the flame propagation mode 
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Fig. 1. Flame propagation into quiescent fuel of a density of 5 x 10 7 gem -3 . Each 
contour marks a time step of 3.3 x 10 -3 s. The contours are shifted for better 
visibility and the spacing between them does not reflect the flame velocity. 

in the context of the delayed detonation model of SNe la. A detonation wave 
is mediated by shocks and propagates with supersonic velocities. 

3 Results of numerical simulations 

In order to find answers to the above questions, we performed numerical 
simulations of the flame propagation in two dimensions. For a description 
and tests of the numerical implementation we refer the reader to [10, 13]. 
Our simulations addressed the flame propagation into quiescent fuel as well as 
flame interaction with a vortical flow field. The simulations were performed on 
an equidistant cartesian grid. The flame was initialized in the computational 
domain and perturbed in a sinusoidal way from its planar shape. To follow 
the long-term flame evolution, we performed the simulations in a frame of 
reference comoving with the flame. This was implemented by imposing an 
inflow boundary condition ahead of the flame front and an outflow boundary 
on the opposite side of the domain. Transverse to the direction of flame 
propagation we applied periodic boundary conditions. 

3.1 Flame propagation into quiescent fuel 

The first part of our investigations addressed the flame propagation into qui- 
escent fuel. Here two stages of flame evolution are revealed. In the beginning, 
the small sinusoidal perturbation imposed on the flame front increases due 
to the LD instability. This part of flame evolution was studied in [10], were it 
was shown that the growth rates observed in our simulations agree well with 
the expectations from Landau's linear stability analysis [6]. In the second, 
nonlinear stage the flame stabilizes in a cellular pattern. The propagation of 
a flame with an initial perturbation of 6 periods fitting into the domain is 
plotted in Fig. 1. The tendency of the small cells to merge in the nonlinear 
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Fig. 2. Flame propagation into vortical fuel of a density of 5 x 10 7 gcm 3 . (a) 
v' / si = 0.7; (b) v' / si = 2.5. Each contour represents a time step of 8.0 x 10 -4 s (a), 
and 2.4 x 10~ 3 s (b). A gain, the contours are artificially shifted for better visibility. 

regime finally stabilizing in a single domain-filling cusp-like structure is ap- 
parent here. This result is consistent with semi-analytical studies of flame 
evolution [14]. 

3.2 Flame propagation into vortical fuel 

In SN la explosions, turbulent velocity fluctuations around the Gibson scale 
can be expected and also relics from pre-ignition convection may contribute 
to turbulent flows. In order to explore the effects of flame interaction with 
such velocity fluctuations, we modified the inflow boundary condition to a 
so-called oscillating inflow condition, which generates vortices by imposing 
the following velocities at the boundary 

v x = si { — 1 + v sin 2kTry cos 2kn{x — tsi ) } 
v y = s i v' cos 2k-ny sin 2kn(x — tsi ) . 

Here, the parameter v' characterizes the strength of the imprinted velocity 
fluctuations and k denotes the wavenumber of the oscillation. This produces 
what is termed "square vortices" by [15]. 

The results of two of such simulations for different strengths of the square 
vortices are plotted in Fig. 2. Although in case of weak vortices in the in- 
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coming flow there is some interaction with the flame, the flame still shows 
the tendency to align in a large- wavelength cellular structure (cf. Fig. 2a). 
This is similar to the case of flame propagation into quiescent fuel. On the 
other hand, if the vortices are strong enough, they can completely break up 
the cellular stabilization and prevent the flame from evolving into a single 
domain-filling structure (cf. Fig. 2b). However, we do not observe a dras- 
tic increase in flame surface and thus flame propagation speed. The flame 
structure rather adapts to the vortices imprinted on the fuel flow. 

4 Conclusions 

The presented simulations of flame propagation into quiescent fuel are in good 
agreement with theoretical expectations. The linear stage of flame evolution 
is consistent with Landau's dispersion relation [10, 6] and in the nonlinear 
regime the flame stabilizes in a cellular pattern. Thus, our hydrodynamical 
model of flame evolution shows that the cellular burning regime exists for SN 
la explosions. 

Interaction of the flame with vortical flow fields may lead to a break-down 
of cellular stabilization if the velocity fluctuations are sufficiently large. In this 
case, however, the flame shows the tendency to adapt to the imprinted flow 
field. Thus, we observe a moderate increase in the effective flame propagation 
velocity but no drastic effects. No indication of active turbulent combustion 
could be found. Therefore, it seems unlikely that effects around the Gibson 
scale account for a DDT at late stages of the SN la explosion as has been 
anticipated by [11]. More detailed discussions of flame propagation into qui- 
escent fuel at different fuel densities and interaction with vortical flow fields 
of varying strengths are currently in preparation [16, 17]. 
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